Fractal structure of near-threshold quarkonium production off cold nuclear matter 
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We investigate near-threshold production of quarkonium resonances in cold nuclear matter 
through a scaling theory with two exponents which are fixed by existing data on near-threshold 
J/V' production in proton-nucleus collisions. Interestingly, it seems possible to extend one of the 
multifractal dimensions to the production of other mesons in cold nuclear matter. The scaling theory 
can be tested and refined in experiments at the upcoming high-intensity FAIR accelerator complex 
in GSI. 
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There are unexplored systcmatics for the production of 
quarkonia close to threshold. For the J/ip, a few exper- 
imental studies were carried out in the days before the 
understanding of QCD was mature. Attention shifted to 
the high-energy frontier, since perturbativc QCD turned 
out to be the tool appropriate to that region, and weak- 
coupling theory was applied successfully to higher-energy 
production of quarkonium in pA collisions [1HJ|. How- 
ever, the planned FAIR in GSI presents a grand opportu- 
nity for the study of near-threshold particle production in 
cold nuclear matter, and test and refine a scaling theory 
which is developed here. This is interesting for several 
reasons. First, as we argue here, testing the limits of 
such a scaling theory shows where the crossover between 
hadron and quark descriptions of matter lie. Second, this 
measurement can be used to test factorization where it 
cannot be proved, and therefore has implications on the 
understanding of CP violations. Third, a detailed under- 
standing of cold nuclear effects in quarkonium production 
is important to tests of the formation of the QCD plasma 
in heavy-ion collisions. Finally, the discovery of scaling 
laws, which is our main result, is of fundamental interest 
since the exponents can define universality classes across 
very broad ranges of physical phenomena. 

The basic toolkit which we bring to this study is the 
modern understanding that the renormalization group 
is expressed via scaling. We wish to ascertain whether 
there is a dynamical symmetry of cold nuclear matter 
so that a small change in the amount of nuclear matter 
can be compensated for by a corresponding change in 
the energy of the probe. Invariance under such scaling 
would manifest itself in the form of exponents which are 
eigenvalues of renormalization group transformations [5[ . 
These scaling exponents are also called fractal dimensions 
or anomalous dimensions. At high energies they can be 
computed in perturbation theory. Near the threshold of 
particle production, they have to be discovered in data 
and understood non-perturbatively. Discovery of scaling 
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exponents, which is our main result, conversely implies 
the scaling symmetries. 

We are interested in the production of a quarkonium 
state, H, in pA collisions near the threshold energy v^oi 
which is the minimum energy in the center of mass re- 
quired to produce H. We follow the convention of writing 
the CM energy, vo, in the equivalent pp system; for a 
fixed target configuration this means S = 2Mp{Ei, + Mp), 
where Et is the beam energy and Alp the proton mass. 
The threshold energy, \/S() ~ 2Mp + M„, where M„ is 
the mass of H. The total inclusive cross section, a, can be 
a function of VS, M„, Alp, and the nuclear mass, AIa- 
Then, a dimensional argument allows us to write 



5oCT = f{A,Y,h), where A 

'-''Mi'-" 



^ Ma 
Mp' 
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and / is a dimensionless function. In this definition of 
A, we neglect the effects of nuclear binding, which are 
expected to be less than 1%, and isospin effects, which 
could be slightly larger. We take masses and branching 
ratios from [6|. 

In this paper we report a scaling analysis of J/ip cross 
sections in a dilepton channel in pp and pA collisions 
from the lowest up to ISR energies J7H22|. but not be- 
yond. Within this data corpus, corrections for kine- 
matic acceptance limitations of each experiment needed 
for global analyses are discussed in [3|, |23[ . We decided 
to examine Ba rather than a, where B is the branching 
ratio in the diclcctron or dimuon channel. The reason 
is that over the years the value of B has moved by more 
than its error bar. When the inclusive cross section in one 
of these dilepton channels is measured, this uncertainty 
does not affect the result. Some experiments correct their 
data for nuclear effects according to a formula A", with 
a obtained from their data. We undid this correction, 
since this is part of our global analysis. 



Near threshold the variable Y is small and close to 
zero. In proton nucleus scattering, we expect some Fermi 
motion: even though the CM of the nucleus is at rest, 
individual nucleons may be moving. The typical energy 
of this movement is of the order of the binding energy 
per nucleon [2J|, and hence comparable to other effects 
which we have neglected. Clearly, Fermi motion can be 
detected with experiments close to threshold since the 
cross section vanishes otherwise. However, for Y > 0.1, 
the effect can be neglected. 

When -s/S" is large, then one expects to be able to com- 
pute f{A, y, h) in models inspired by perturbative QCD; 
this is true whether h is large and Y small [25| or h is 
small and Y large \u\M- ^^ important pre- requisite for 
these computations is the factorization of the initial state 
into parton distribution functions. If these factorization 
theorems were valid then certain kinematic scalings could 
be expected to hold 26| which are seen to fail for F w 1 
22l |. So, in the near-threshold region for J/^j perturba- 
tive QCD inspired models for quarkonium production is 
not viable. 

When Y is small, one expects the cross section to have 
a Taylor expansion in Y. In pp collisions, where A = 1, 
one may then expect the particularly simple parametriza- 
tion fi{Y,h) = Y^fi{h), where /? is a scaling dimension 
and h is a scaled variable. Such a scaling law indeed 
describes the A = 1 part of the data corpus well with 

S/i/S'o = 2.0±0.4nb, /? = 3.20 ± 0.26, (2) 

with covariance of —0.934. 

More generally, we can write f{A,Y,h) = Y^f{A,h), 
where the "renormalized" variables are taken to have the 
form A = A/{y+Y)^, where fj, is an anomalous dimension 
and y is an additive renormalization constant. Then a 
change from Y to Y' = (Y can be compensated by a 
change of the target nucleus from ^ to ^' = S^A with 
^1/m _ (^(^Y + y) / (Y + y) . One has a similar renormalized 
h{h,Y). A power law f{A,h) = A^f(h) then gives a 
multifractal exponent ^"^^^ . Generally, such a power law 
parametrization may be expected when A ^ 1. Using 
the extensive data taken by NA50 [20| we find that this 
is indeed a very good description of the data for A > 
50. Furthermore it is significantly less probable that this 
behaviour extends to all A. 

Moreover, from the data corpus it is clear that a con- 
stant power, independent of Y , is inadequate. Instead 
we write. 



fiA,Y,h) = Y^A''^^^fih). 



(3) 



and choose a particularly simple multifractal exponent 
with a linear dependence on Y. The data corpus supports 
such a scaling law with 
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FIG. 1: The multifractal exponent a for the production of 
J/^p in pA collisions. The band encloses the 68% confidence 
limits of the two parameter fit of eq. Q. 



with covariance —0.984, as shown in Figure [TJ The data 
from CERN-PS [lO| was not used in this fit, because of 
the large errors in this measurement; its inclusion does 
not change the central values of the fit significantly but 
increases the errors. It would be interesting to test in fu- 
ture experiments whether this exponent is modified near 
y = 0.1 as the thresholds for V'(2S) and the xs are ap- 
proached. 

For a test of the scaling behaviour in eq. ([3]) we con- 
structed the scaled cross section Ba/A" , using the largest 
A from each experiment which had A > 50. Fitting by a 
power of Y gives 



B/yS-o = 3.2 ± 0.5 nb, /3 = 3.0 ±0.3, 



(5) 



a = (0.76 ± 0.02) + (0.10 ± 0.01)r, 



(4) 



with covariance of —0.970. In dimensionless units we 
have / = (3.4 ± 0.5) x 10"'^. The exponent is consistent 
with that obtained in eq. (O. In fact, we have plotted 
the pp data also in Figure [2] to show that the two sets of 
data are close, but not identical. This is consistent with 
our expectation that the multifractal exponent of A is 
valid for large A. It is clear from Figure [2] that the data 
corpus does not provide a very stringent test of scaling 
in near-threshold production cross sections. Additional 
data would therefore be very useful. 

It is interesting to extend the parametrization of eqs. 
PI4I5P to other resonances. The exponent a for the 
ground state quarkonium is unlikely to have a close re- 
lationship to that for higher resonances, since they are 
more susceptible to decay due to medium effects. On 
the other hand, one may ask whether there is any uni- 
versality in exponents of ground state quarkonia. To an- 
swer this, we examine the T, for which E772 [27[ re- 
ported a = 0.962 ± 0.006 ± 0.008 at v^ = 38.8 GcV. 
Extrapolating the fit in cq. ^ gives a == 0.883 ± 0.008. 



To improve this we add a resonance mass dependence: 
a = ttQ + aiY + ahh. The fit of eq. (j4]) along with 
the above measurement gives a^ = 0.012 ± 0.002, and 
uq = 0.64±0.02. Interestingly optical models of shadow- 
ing in very low-energy nucleon-nucleus scattering predict 
ao = 2/3. This value is in agreement with results for low 
energy tt production [28|. Coincidentally, the formula is 
also in agreement with the measured value of a for K, p, 
and w production at low energy [29|. 

It is interesting to compare the parametrization of eq. 
([3]) with other studies. The Glauber model treated in the 
eikonal approximation has also been used in cold nuclear 
matter with the parametrization 



f{A,Y,h)=Aexp(^-jA'/') f{Y,h). 



(6) 



In the model, the dimensionless number 7 = paats^, with 
p being the nuclear density, aats having the interpreta- 
tion of a cross section for absorption of J/ip in cold nu- 
clear medium, and XA^^"^ being the path length of the 
J/ip in the nucleus (see, for example, [30|). Although the 
functional form F{A) seems to be very different from A", 
they can be numerically close. As a result, the Glauber 
model cannot be experimentally distinguished from the 
simpler model of eq. ([3]). In fact, the data corpus can be 
used to parametrize 7 = (0.27 ± 0.02) - (0.092 ±0.007)y. 
Then, using the values p = 0.16/fm^ and A = 1.1 fm, we 
can write 



(^abs 



(15 ± 1) - (5.2 ± OA)Y mb, for Y < 2.5. (7) 
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This is consistent with the values extracted in 

Different forms have also been used for the scaling 
function .f(Y , h) = SoaA~". A parametrization sug- 
gested in [si! is compatible with the data for F > 1, but 
differs from that of eq. ([3|) near threshold. A parametriza- 
tion of the threshold effect as (1 — 1/ expl")'^. as in [16| . 
is compatible with data for Y > 1, but differs from eq. 
((HI) near the threshold, as shown in Figured Fermi mo- 
tion based explanations of the discrepancy with data for 
Y < 1 are ruled out because two pieces of the data come 
from pp experiments. 

We argue that the crossover between these forms is 
physical and interesting. Near threshold y ~ (1 ~ 
1/expy), so the powers /3 and v can be compared. 
The fit value /3 ~ 3 can be interpreted within the old 
"spectator counting rules" [32], as indicating that there 
are 2 spectators. Since hidden charm production is an 
OZI-violating process, this is consistent with a purely 
hadronic origin of near threshold processes. On the other 
hand, the best fit value of z/ ~ 11 (Figure [2]) implies that 
there are 6 spectators. This is consistent with counting 
the valence quarks from two participating hadrons in a 
process dominated by gluon fusion. So the crossover from 
one regime to another may give us a clue to the regime 
of validity of perturbative QCD. 
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FIG. 2: Scaling plot for the experimental data on the total 
inclusive cross section for the production of J/tj) in pA colli- 
sions. The band enclosed by full lines is the 68% confidence 
limit of the two parameter fit of eq. Q to pA data for A > 50, 
so excluding p-Be data from IHEP @] and E331 [y| and p-Si 
data from E771 [10]. The band enclosed between the dashed 
curves is the 68% confidence limit of a two parameter fit to 
the form K[l — e"^]" to all pA data. 



The scaling theory is unable to relate the values of f{h) 
in experiments with different initial states: for example 
pA, pA, TT^ A, etc. If some form of the factorization the- 
orems were valid, then there could be relations such as 



fjpA -^H) ^ /{ttA ^ H) 
fivA'^H) fi^A'^HY 



(8) 



where f{ab — > h) — uSq and a is the total inclusive 
scattering cross section for the reaction a6 — J- /i at a fixed 
Y . Such tests could also be performed with light ions 
instead of pions. The validity of factorization would be 
an interesting investigation especially in view of the fact 
that the connection between hadron decays and the CKM 
matrix elements can only be made with this assumption 
near threshold. 

Proceeding beyond total inclusive cross sections near 
threshold is hard due to the paucity of data. There is 
scattered evidence for cold nuclear effects in (pj), from 
NA38 [11] and HERA-B [2lj. However, there is too lit- 
tle data for a systematic study of the effect. The sparse 
corpus shows a roughly linear rise of (p^) with Y , from 
a vanishing value at y = 0. Clearly, high statistics stud- 
ies of Pt and Xp distributions of the J /ip near threshold 
would be very welcome. 

Since data with beam energy E}, < 100 GeV is very 
sparse [3,[l0:ll4|i ^^^ SIS-100 accelerator at GSI presents 
an opportunity to probe the region of y < 0.4 very thor- 
oughly with modern statistics. With a beam luminosity 
of 1 Hz/nb; fair event rates could be obtained, and these 



scaling laws can be tested well. This would make the GSI 
an ideal test bed for exploring the near-threshold pro- 
duction process for J/ij) as well as cold nuclear effects, 
including questions about factorization and p^, and Xp 
distributions. A wish list would contain measurements 
with pp and a variety of pA collisions to check the scal- 
ing of eq. dl]). The pp data could also be used to check 
the scaling exponent of eq. ([2]) and whether it is compat- 
ible with the pA result of eq. ([5]). A range of A can be 
used to test the region of validity of the power law A"" . 
A systematic study of p^. and Xp distributions would also 
be extremely useful. 

In summary, we have extracted a power law 
parametrization of the cross section for near-threshold 
J /ij) production off cold nuclear matter. The results are 
given in eqs. pi4l5p . Such power laws are more than just 
a parametrization, since they reveal certain dynamical 
symmetries of hadronic systems, which equate a physi- 
cal system with one Y and A to another with different Y 
and A. These identities constitute renormalization group 
transformations, and should eventually be computable 
from QCD. Interestingly, it seems that the multifractal 
exponent a of eq. (|4]) can be extended to the production 
of the T via 

a{Y,h) = (0.64 ±0.02) + (0.10 ±0.01)r 

-^-(0.012±0.002)/i. (9) 

Coincidentally, this form it also reproduces the exponent 
required for inclusive tt, K^ p and uj production, but not 
for ip' or (f) production. These scaling laws present fun- 
damental restrictions on QCD, and therefore should be 
of priority in upcoming low-energy and high-intensity ex- 
periments at the SIS-100/300. 
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